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INTRODUCTION 


This  analysis  will  attempt  to  predict  propellent  pressures  which 
act  on  various  parts  of  a muzzle  brake.  These  pressures  will  then  serve 
as  inputs  to  a stress  analysis  of  the  brake.  The  105mm  M68  gun  and  an 
experimental  howitzer  brake,  fitted  to  the  gun,  were  selected  because 
the  gun  is  being  fired  in  a life  test.  For  economy,  pressures  on  the 
brake  will  be  measured  during  a few  of  those  rounds  to  check  the  analysis. 

The  internal  pressures  created  in  approximately  the  first  66  micro- 
seconds following  the  barrel  uncorking  are  analyzed.  The  analysis  is 
broken  into  two  parts,  one  for  the  diffuser  section  of  the  brake,  and 
one  for  the  downstream  portion  through  which  the  gases  exhaust  to  the 
atmosphere.  Both  flows  are  modeled  as  one-dimensional  during  the  short 
period  of  time  which  it  takes  for  the  uncorked  gases  to  reach  the  end 
of  the  brake. 

No  attempt  is  made  to  predict  the  flow  fields  or  the  pressures 
produced  by  it  at  subsequent  times,  although  approximate  methods  might 
be  developed  to  do  so,  should  the  predictions  of  this  analysis  yield 
results  in  reasonable  agreement  with  experiment. 

ANALYSIS  OF  BLOW-BY 

At  the  instant  that  the  projectile  uncorks  the  barrel  a shock  wave 
is  formed  in  the  brake  at  the  muzzle.  This  shock  wave  propagates  down 
the  diverging  annular  space  between  the  projectile  and  the  inside  wall 
of  the  muzzle  brake.  The  flow  initiated  by  this  shock  wave  represents 
a leakage  of  high  pressure  gas  from  behind  the  projectile.  (See 
Figure  1 ) . 
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Figure  1.  Projectile  and  brake  at  muzzle  uncorking. 


Figure  2 shows  the  various  regions  associated  with  the  sudden 
release  of  high  pressure  gas  into  the  annulus  following  the  uncorking. 
The  gas  which  escapes  through  this  annulus  traverses  an  area  change 
with  a ratio 

0)  At  = (5)2-(4)2  = 9 

A4  (4)2  16 

The  equivalent  diameter  ratio  for  this  area  change  is 

(2)  D]  = \ = 3 

D4  \j  A4  4 

According  to  Alpher  and  White1  the  effect  of  an  area  change  in  a 
shock  tube  can  be  accounted  for  by  introducing  a factor  (g)  into  the 
standard  shock  tube  "burst"  equation  so  that  it  takes  the  modified  form 

(3)  1 f 2"^  Ms2  - (Tj-n' 

p4  = g L 7\  + l _ 

Pl  >4-1  2^4 


In  this  equation  subscripts  refer  to  the  regions  so  numbered  in  Fig  2 
and: 

Tl  = 530°R 

T4  = 31 50°R 

Y]  =1.4  (air) 

Y *4  - 1.25  (combustion  products) 

1.  R.  A.  Alpher  and  D.  R.  White,  "Ideal  Theory  of  Shock  Tubes  with 
Area  Change  Near  the  Diaphragm"  G.  E.  Rpt  57-RL-1664,  Jan  1957 
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V4  = 3800  ft/sec  (projectile  velocity) 
g = 1.25  for  D4/D]=  4/3 

Rl  = 53.3  ft-1 bf/lbm  °R  (air) 

R4  =67.0  ft-1 bf/lbm  °R  (combustion  products) 

P 4 = 6.82  lbm/ft3 


gc  = gravitational  constant  = 32.2  Ibm  ft/lbf  sec^ 
C4  =^4"  9c  T4  = 2918  ft/sec 
P4/P1  = 10,000/15  = 667 


Cl 

c4 


l/l.4  (53.3)  530 
7 1.25  (67)  3150 


3.86 


With  these  inputs,  equation  (3)  becomes 


(3a) 

667  = 0.933MS2  - 0.133 

[l.  163-0. 0393(MS-1/MS)]10 

Solving  this  equation  (3a)  for  the  initial  shock  Mach  number  (Ms) 
yields  a value  of  9.2.  The  shock  wave  itself  travels  with  a velocity 
of  9.2  (1130)  = 10,396  ft/sec  while  the  fluid  which  it  sets  into  motion 
moves  with  a velocity  given  by 


(4) 


A (1130)  (9.2  - 1 ) 
2,4  ( 9.2) 


V2  = 85.61  ft/sec 

Obviously  the  shock  wave  and  the  fluid  behind  are  traveling  much 
faster  than  the  projectile, reaching  the  end  of  the  diffuser  section  of 
the  brake  (4"  from  the  muzzle)  before  the  base  of  the  projectile  does. 
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The  propagation  of  a shock  wave  through  an  area  change  has  been 


investigated  by  Whitman2  who  showed  that 

Koo 


(5) 


Me  = Ms  (As) 

(A?) 


where  .36 


^ 9.2 


(5)2-  (4)2 
(7)2-  (4)2 


0.18 


Me  = 7.28 

The  conditions  behind  the  shock  wave  at  the  exit  of  the  diffuser 
are  given  by 

(6). 


!i  = 2?1  2 $C± 

Pi  ^1+1  Me  ‘^+1 


(7) 


(8) 


YJ  (7.28)2 


0.4 

2.4 


Pe 

Pi 


p7  = 61.7 


r 

r 

Te  - 

Tp  = 


(fi+1)  MP2 

(/l-l)  Mg2  + 2 


5.48 

Pe  P\ 

PTTe 

§7^t  530  = 5964°R 


Tl 


2.  G.  B.  Whitman,  "On  the  Propagation  of  Shock  Waves  through  Regions 
of  Non-Uniform  Area  Flow,"  J.  Fluid  Mechanics,  Vol . 2,  P.  337  (1958) 
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(9) 


v 


e 


ve  = 6726  ft/sec 


The  leakage  passed  by  the  projectile  (blow-by)  can  be  estimated  from 


The  time  for  the  projectile  base  to  reach  the  end  of  the  diffuser  is 

- = 4/12  = 8.77  (10-5)  sec  so  that  the  total  leakage  during  this  period 

x 3800 

is  0.0436  Ibm  which  is  only  a few  percent  of  the  total  mass  and  therefore 
quite  negligible. 

ANALYSIS  OF  THE  GAS  FLOW  BEHIND  THE  PROJECTILE 

Since  very  little  of  the  barrel  gas  is  lost  via  blow-by  as  the  pro- 
jectile base  moves  through  the  diffuser  section,  most  of  the  gas  which 
was  behind  the  projectile  base  at  uncorking  must  still  be  behind  the 
base  when  it  reaches  the  diffuser  exit  (e). 

If  it  is  assumed  that  this  gas  undergoes  an  isentropic  expansion, 
then  the  pressure  at  this  point  must  be  approximately 


(10) 


144  4 

m = 497.6  lbm/sec 


m = A ve 

= 5.48(0.075)  (49-16)  (6726)77 


Therefore  pJ  = 10,000  (1_6) 


(49) 

P41  = 2468  psia 

02) 
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03) 


(14) 


C,  = 


4-  4 


4l 

C 1 = Jl. 25(32. 2)67(2381) 
4 

C.i  = 2534  ft/sec 
4 


f4i  ■ ^4  > 


(‘V  - 6. 


V41 

82  (16)  = 2.23  lbm/  „ 
(49)  ft^ 


ANALYSIS  OF  DIFFUSER  UNCORKING 


Figure  3 illustrates  the  situation  assumed  to  exist  at  the  instant 
that  the  projectile  base  reaches  the  end  of  the  diffuser.  It  will  be 
assumed  that  the  disk  of  fluid  shown  at  the  left  of  this  figure  propa- 
gates to  the  right  with  the  projectile  velocity  x,  (3800  ft/sec)  while 
expanding  radially  behind  a radial  shock  wave.  This  assumption  is 
analogous  to  that  used  in  the  report  by  Soifer^  with  the  exception  that 
the  explosion  does  not  emanate  from  a line  source  and  is  one-dimensional 
rather  than  cylindrical.  The  assumption  of  one-dimensional  flow  seems 
more  realistic  considering  the  presence  of  the  brake  side  plates  beyond 
the  diffuser  section.  The  burst  of  a finite  size  volume  assumed  here 
avoids  the  rather  unrealistically  low  densities  in  the  core  of  the  burst 
which  result  from  the  assumption  of  a line  source. 

The  shock  Mach  number  for  a shock  wave  propagating  radially  from  the 
"window"  of  the  muzzle  brake  is  given  by  the  equation 


3.  "Muzzle  Brake  Analysis,"  S & D Dynamics,  Inc,  TR  # 73-6,  Oct  73, 
Contr.  No.  DAAF  07-73-0541 
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05) 


& _ 


mA djzhuA 

YL±JL 


(15a) 


2468 

/S' 


//-,04t4fi7s-A 


j ~j/o  / <£  ■4- 

Ts)J 


This  equation  is  the  same  as  equation  (3)  without  the  factor  (g). 
Furthermore,  the  last  term  in  the  denominator  of  (3)  vanishes  here  because 
v4l » radial  velocity,  is  initially  zero  in  this  radially  expanding 

shock. 

The  shock  Mach  number  for  this  one-dimensional  shock  wave  moving  at 
a right  angle  to  the  flight  path  is  Ms  = 4.25.  Behind  this  shock  wave 


(16) 

07) 

(18) 

(19) 


u2'  = Jrf+T  Cl 


(M 


's  M, 


) ~ 3774  ft/sec 


,jl=  «>Vl)Ms2 

p (rrl)  Ms2  +2  = 5.40 

P21  - ztf  Ms2  -(Tr1 
STfi- 


= 20.9 


Tol 


f2l.fl 

’l  fil 


T-|  = 2051  °R 


The  time  required  for  the  projectile  base  to  move  from  the  end  of  the 


diffuser  to  the  baffle  is  x _ 3/12  6.58  (10-5)  sec. 

* 3800 

During  this  time  the  leading  edge  of  the  radial  shock  wave  travels 
a distance 


(20)  A R = C-jMs  x 

x 

A R = (C-] Ms ) 3 = 3>78  in 
(3800) 
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as  shown  in  Figure  3.  The  trailing  edge  of  the  centered  expansion  fan 
travels  a distance 


(21) 


A r = 04!  x.  = 04!  t 


x 


= 2.00  in 


Figure  4 shows  the  relative  position  of  the  characteristics  and  the 
shock  front  with  respect  to  the  baffle  when  the  moving  disk  of  gas 
impinges  on  it. 

In  order  to  find  the  pressure  distribution  on  the  baffle,  the 
"backward"  running  characteristics  (shown  as  a fan  in  Figure  4)  are  used 
to  establish  the  location  on  the  baffle  while  the  forward  running  char- 
acteristics are  used  to  carry  the  core  conditions  through  the  fan. 

For  example,  the  characteristic  which  intercepts  the  outer  edge  of 

the  baffle  has  a slope  of  1.5/12  _ _ n0nn 

'6758 (TO-5)  " 1900  tt/SeC 

therefore,  the  slope  is  equal  to  U0-C0  = 1900  where  the  subscript  0 refers 
to  the  outer  edge  of  the  end  plate.  Since  along  the  "forward"  running 
characteristic  (not  shown  in  Figure  4),  the  value  U-/-2  c = constant  = 


8C4I , i.e. , 


f 4-1 


(23)  U0+8  C0  = 8 C41  = 8 (2534) 


Combining  these  two  equations  yields 


C0  = 2041  ft/sec 
U0  = 3941  ft/sec 


and 


from  (13)  7“  - QclL 


o 3/,zrR 
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Figure  4.  Wave  diagram  for  shock  wave  propagation  through  the  window. 


£83  fsia. 


(24) 


The  stagnation  pressure  rise  due  to  momentum  exchange  as  the  disk 
impinges  on  the  baffle  is 


The  total  pressure  felt  at  the  extremity  of  the  baffle  is  therefore 


Using  a similar  procedure  the  distribution  of  total  pressure  can  be 
mapped  for  the  entire  baffle.  The  following  sections  will  describe  the 
results  of  calculations  for  the  pressure  distribution  throughout  the 
brake  interior. 

INITIAL  PRESSURE  DISTRIBUTION  IN  THE  4"  LONG  DIFFUSER  SECTION 

The  pressure  distribution  imposed  radially  upon  the  diverging  dif- 
fuser section  by  the  isentropic  gas  expansion  directly  behind  the  base 
of  the  projectile  is  computed  from  (11)  and  summarized  in  the  following 
table. 


(25) 


(26) 


P0  total  “ 1227  + 283  1510  psi 


Diffuser  Diameter  (in) 


Table  1 


Radial  Pressure  (Psi) 


Following  the  "uncorking"  of  the  diffuser  by  the  projectile  base, 

the  diffuser  pressures  will  decrease  rapidly  tending  toward  the  steady 

values.  As  a consequence  the  values  given  in  Table  1 should  represent 

the  maximum  pressures  experienced  by  the  diffuser. 

INITIAL  PRESSURE  DISTRIBUTION  ON  THE  SIDE  PLATES  BETWEEN  THE  DIFFUSER 
AND  BAFFLE 

Following  the  uncorking  of  the  diffuser  the  results  of  a one- 
dimensional shock  expansion  from  the  "windows"  of  the  brake  yield  the 
following  initial  side  plate  pressure  distributions: 

Table  2 

(1"  from  Diffuser  Exit) 


Distance  from  Center  Line  (in)  Pressure  (Psi) 

0 - 2.83"  2468 

3.0  1864 

3.5  759 

4.0  284 

4.5  95 

5.0  28 

Table  3 

(2"  from  Diffuser  Exit) 

Distance  from  Center  Line  (in)  Pressure  (Psi) 

0 - 2.17  2468 

2.5  I860 

3.0  1202 

3.5  759 

4.0  470 

4.5  285 

5.0  167 
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Table  4 

(3M  from  Diffuser  Exit) 

Distance  from  Center  Line  (in)  Pressure  (Psi) 


0 - 1.5 

2468 

2.0 

1860 

2.5 

1391 

3.0 

1035 

3.5 

759 

4.0 

552 

4.5 

399 

5.0 

283 

The  location  of  the  pressure  stations  cited  in  Tables  2-4  are  shown 
in  Figure  5. 

INITIAL  PRESSURE  DISTRIBUTION  ON  THE  BAFFLE 

The  total  pressures  which  are  imposed  initially  by  the  arrival  of 
the  gases  directly  behind  the  projectile  base  are  shown  in  Table  5.  The 
location  of  the  pressure  stations  used  are  shown  in  Figure  6.  These  are 
calculated  by  the  method  shown  previously  using  equations  (22)  through  (26). 


Table  5 

Distance  from  Center  Plane  (in) 

Pressure  (Psi) 

0 - 1.5 

9413 

2.0 

7397 

2.5 

5782 

3.0 

4501 

3.5 

3462 

4.0 

2648 

4.5 

2015 

5.0 

1510 

15 

Pressures  given  at  these  grid  points 


5: 
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Figure  5.  Location  of  side  plate  pressure  stations. 


Pressures  uniform  verticil 1y 
along  these  lines  — — 


Figure  6.  Location  of  baffle  pressure  stations. 


These  values  are  the  stagnation  pressures  obtained  by  summing  both 
the  static  and  dynamic  pressures. 

All  interior  surfaces  of  the  brake  will  experience  some  pressures 
prior  to  those  listed  above  due  to  the  precursor  flow  caused  by  blow-by. 
These  pressures  are  expected,  however,  to  be  masked  in  time  and  amplitude 
by  the  almost  simultaneous  arrival  of  the  uncorking  flows. 

SECOND  STAGE  ANALYSIS 

The  techniques  described  above  can  be  extended  to  the  analysis  of  a 
second  stage  of  the  muzzle  brake  with  some  slight  modifications. 

Figure  7 shows  the  wave  diagram  for  the  radial  motion  of  gas  in  the 
second  stage  as  the  disk  of  gas  expands  behind  the  projectile  base.  The 
only  backward  running  characteristics  generated  in  the  first  stage  which 
pass  into  the  second  stage  are  those  directly  behind  the  projectile  base 
when  it  passes  through  the  first  baffle.  If  it  is  assumed  that  this 
baffle  is  quite  thin,  then  the  burst  of  gases  into  the  second  stage  gives 
rise  to  a new  centered  expansion  fan  as  shown  in  Figure  7. 

From  the  previous  calculations  of  the  values  in  Table  4,  the  condi- 
tions two  inches  from  the  center  line  at  the  moment  that  the  projectile 
passes  through  the  first  baffle  are 

C411  = 2463  ft/sec 
T4»  = 2250  °R 

(27)  f 4.,  = 1.778  lbm/ft3 

P411  = 1860  psia 
U411  = 563.6  ft/sec 

The  shock  Mach  number  for  the  radial  growth  of  the  disk  after  it 
passes  through  the  first  baffle  plate  is  given  by  (28)  which  is  similar 
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Figure  7.  Wave  diagram  for  second  baffle. 


- /./67M?- 0.22  8/2 
p,  02  g/2  - a,  0478(%  ” 


'J 


The  shock  Mach  number  which  satisfies  this  equation  is  Ms  = 4.3. 
The  velocity  of  the  gas  behind  this  shock  wave  is  from  (4) 


(4) 


(29) 


^2"  = _? Ms-1 

Y7+1  C1  ( Ms) 

= f-r  (1130)  (4.3-1  ) 

4,3 

U2 m = 3830  ft/sec 

The  pressure  in  this  region  is  from  (6). 

P2»  = 2&1  Ms2  - (^1-1)  = 21.4 

P1  Fy+ 1 )^+l 


M ~ 321  psi 

If  the  distance  between  the  first  and  second  baffles  is  assumed  to 
be  3.0  inches,  then  the  shock  wave  will  have  advanced  radially  a distance 
equal  to 
(30) 


sCl 


1130(4.3) 


3j(3800j/=  3((  3800  )J=  3.84" 

which  is  0.84"  beyond  the  edge  of  the  second  baffle  at  this  point, 
assuming  that  the  second  baffle  is  the  same  size  as  the  first  baffle. 

The  characteristic  which  forms  the  left  side  of  the  expansion  fan 
coincides  with  characteristic  along  which  U4"  and  04-1  have  constant 
values  of 
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from  (27). 


U4"  = 563.6  ft/sec 
C4"  = 2463  ft/sec 


Regions  4"  and  3"  are  connected  by  a forward  running  characteristic 


along  which 


(31) 


U3" 


+ 2 
T47' I 


C3n 


563.6  + 8(2463)  = 3830  + 8C3«.  ' 


C3"  = 2055  ft/sec 


Then  from  (13) 
(32) 


T3"  = 1566  °R 


from  (12) 
(33) 


(34)  from  (8) 


P3"  = 304  psia 
^3"  = (• 4"  P3"  J4" 

P4"  T3 11 

=1.778304  2250 

1860  1566 


P 3"  = 0.132 

and  from  (25)  the  stagnation  pressure  is 

(35)  (p3»)  = 


= 0.132  (3800) 

32.2  VT27 
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P311  = 411  psi 

1 he  total  pressure  in  the  constant  velocity  region  behind  the  shock 
wave  is  therefore 


(36)  304  + 411  = 715  psia 

The  leading  edge  of  the  expansion  fan  has  a slope 


(37) 


= 3830  - 2055 


= 1775  ft/sec 

The  radial  position  of  the  leading  edge  of  the  fan  at  the  second  baffle  is 

(38)  AR1  =(1775)6. 58(1 0"5) (12)  = 1.40  in 
Referring  to  Figure  7 and  (36)  shows  that  the  outer  1.6  inches  of  the 
second  baffle  plate  has  a constant  total  pressure  of  715  psia  impressed 
upon  it. 

The  backward  running  characteristic  coincident  with  the  hole  in  the 
second  baffle  has  a zero  slope  on  the  wave  diagram  so  that  Uf-Cjs0  or 
'■j-C-j  along  this  characteristic.  The  forward  running  characteristic 
which  connects  the  hole  with  the  backward  running  characteristic  at  the 
edge  of  the  fan  carries  a constant  value  of  U+  2 C so  that 

r4-i 

(39)  2 2 

^i+  (*4~l")  ci  = U^i't  (f^-l ) C^i  11  or 
2 

U.j+  (TrO  Ci  = 563.6  + 8 (2463) 

where  Ui  and  Ci  are  the  velocity  and  sound  speed  at  the  edge  of  the 
hole.  Therefore 


9Ui  = 9Ci  = 563.6  + 8(2463) 
Ui  = Ci  = 2252  ft/sec 
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Then  as  previously 
(40) 


(41) 


(42) 


Ti  = T4m  (C1J 


C411 


= 2250 


T-j  = 1881  °R 


■ ■ ® 

l 

nssi ) 

\2250y 


1860 


P-j  = 760  psia 


. 8 

= 1 .778/(760)1 

[moj 

fi  = 0.869  lbm/ft3 


The  stagnation  pressure  at  the  hole's  edge  is  then 

(43)  (Pi)  -Pi  v2 

s 


go 
= 0.869 


32 


69  {38001 

.2  \ 12  I 


(Pi)  = 2706  psi 
s 

The  total  pressure  at  the  hole  of  the  second  baffle  2.0  inches  from  the 
center  plane  is  2706  + 760  = 3466  psia. 

The  pressure  distribution  throughout  the  fan  can  be  computed  in  the 
same  manner  as  for  the  first  baffle.  Comparing  the  pressures  on  the 
second  baffle  with  those  on  the  first  (see  Table  5)  shows  that  the 
stagnation  pressures  on  the  second  baffle  at  the  hole  and  outside  edge 
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are  roughly  one-half  those  experienced  by  the  first  baffle  plate.  The 
total  force  on  the  second  baffle  plate  would,  however,  be  somewhat  less 
than  one-half  the  force  on  the  first  baffle  due  to  the  fact  that  a large 
portion  of  the  second  baffle's  area  is  exposed  to  a constant  pressure  of 

715  psia. 

The  pressures  on  the  side  plates  of  the  second  stage  are  consider- 
ably lower  than  those  on  the  first  stage,  ranging  from  760  along  the 
center  plane  to  304  just  behind  the  shock  wave. 

Portions  of  the  side  plate  behind  the  first  baffle  will  be  beyond 
the  leading  edge  of  the  shock  wave  in  the  additional  66  microseconds 
covered  by  this  analysis.  Pressurization  of  these  areas  at  future  times 
will  result  from  the  complex  flows  which  are  established  following  the 
collision  of  the  leading  gas  disk  with  the  baffle. 

ESTABLISHMENT  OF  STEADY  FLOW 

The  condition  of  steady  state  flow  will  be  attained  when  roughly 
the  backward  running  characteristic  (which  is  actually  swept  downstream 
in  this  case)  reaches  the  end  of  the  diffuser.  The  time  at  which  this 
occurs  can  be  estimated  by  examining 

(44)  U4  - C4  = 3800-2918  = 882 

which  is  roughly  the  speed  with  which  pressure  signals  are  swept  down- 
stream through  the  diffuser.  Steady  flow  is  attained  in  the  time 
interval 

A "5 

(45)  A t = 4712  = 33.3(10  ) sec 

882 

The  projectile  leaves  the  second  stage  baffle  at 

-5 

(46)  t = 10/12  = 22  (10  ) sec 

3800 
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so  that  very  shortly  after  the  projectile  clears  the  muzzle  brake  the 
diffuser  flow  becomes  steady  (assuming  the  muzzle  conditions  are  quasi- 
steady) . 

Steady  flow  calculations  through  an  area  ratio  49/16,  i.e.,  from 
the  entrance  to  exit  of  the  diffuser  section  show  that  the  steady  state 
exit  conditions  are  as  listed  in  Table  6 along  with  those  for  Region 
4'  immediately  behind  the  projectile  as  it  uncorks  the  diffuser.  These 
latter  are  given  by  equations  (11)  through  (13). 

Table  6 

Steady  State  Behind  Projectile  Base 


Te 

2068  °R 

2381  °R 

ce 

2364  fps 

2534  fps 

ve 

6146  fps 

3800  fps 

Pe 

1220  psia 

2468  psia 

The  major  changes  which  occur  during  the  transition  period  to 
steady  state  occur  in  the  exit  velocity  Ve  and  the  exit  pressure.  As 
far  as  the  side  plates  are  concerned,  this  change  is  advantageous.  Due 
to  the  increase  in  exit  velocity,  it  is  difficult  to  say  whether  the 
steady  state  is  better  or  worse  than  the  transient  effects  analyzed 
above. 
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